Directional tunneling escape from nearly spherical optical resonators 
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We report the surprising observation of directional tunneling escape from nearly spherical fused 
silica optical resonators, in which most of the phase space is filled with non-chaotic regular trajec- 
tories. Experimental and theoretical studies of the dependence of the far-field emission pattern on 
both the degree of deformation and the excitation condition show that non-perturbative phase space 
structures in the internal ray dynamics profoundly affect tunneling leakage of the whispering gallery 
modes. 



Electromagnetic fields in a uniform dielectric sphere 
can be calculated in much the same way as quan- 
tum mechanical wave functions in a spherically sym- 
metric potential. The sphere exhibits whispering-gallery 
(WG) modes, which are long-lived resonances with elec- 
tromagnetic fields confined near the surface 0. For 
small deviations from the spherical shape, e = {Rmax — 

Rmin) I i^Rmax ~t~ Rmin) 1 (where Rmax^ Rmin <^^^ ^r-? 

maximal and minimal radii) , perturbative treatments |2| 
are routinely employed to infer deformation parameters 
from the splitting of azimuthal degeneracy of the WG 
modes Certain strongly non-spherical resonators, on 
the other hand, can be analyzed with methods from quan- 
tum chaos such as random-matrix theory or periodic- 
orbit expansions 0. Many generic resonator shapes, 
however, fall into a transition regime in which none of 
these known methods apply globally. When entering this 
regime from the perturbative side, calculations may en- 
counter singularities and undefined limits 0. Experi- 
mental studies in this regime can thus provide insight 
into how nature resolves the competition between per- 
turbative and non-perturbative physics, here with the 
resonator shape as a control parameter. 

In this paper we present studies of far-field emission 
patterns and resonance lifetimes of deformed fused-silica 
microspheres. The experimental and theoretical studies 
of the dependence of the far-field emission pattern on 
both the degree of deformation and the excitation con- 
dition show that highly directional tunneling escape can 
occur in microspheres with small deformation (e « 1%) 
and with a large size parameter (kR = 2ttR/X « 785, 
where A « 800nm, R « lOO/im). These results are com- 
pletely unexpected in the ray optics limit or in the earlier 
perturbative wave treatment with e as a small parame- 
ter. The observed emission pattern shows that on one 
hand, the ray model breaks down in predicting the es- 
cape mechanism of WG modes for resonators with small 
deformation. On the other hand, the nonperturbative 
phase space structures predicted by the ray model can 
profoundly affect tunneling leakage of the WG modes. 
Even at extremely small deformation, nonperturbative 
phase space structures still persist and are directly re- 
sponsible for directional tunneling escape. This intri- 
cate interplay between ray and wave dynamics provides 



essential physical insights into properties of weakly de- 
formed WG resonators, especially, for those in which the 
phase space is filled with non-chaotic regular trajecto- 
ries. Directional emission patterns have previously been 
observed in more strongly deformed resonators where a 
significant fraction of the internal rays shows chaotic mo- 
tion 0. 

Deformed microspheres were formed by melting to- 
gether two spheres of similar sizes. The individual 
spheres were fabricated by melting an optical fiber tip 
with a focused CO2 laser beam. The spheres were 
brought into contact and heated until surface tension pro- 
duced a completely convex surface. By carefully control- 
ling the temperature of the glass using the CO2 laser, it 
was possible to repeatedly reduce the degree of deforma- 
tion. For reference, we define the elongated axis of the re- 
sulting prolate spheroid, which also connects the centers 
of the two original spheres, to be the x-axis. The z-axis 
is defined by the remaining fiber stem, which breaks the 
rotational symmetry about the x-axis, making the de- 
formed microsphere completely non-axisymmetric. The 
stem was held by a fiber chuck for easy manipulation of 
the sphere. Images of these deformed microspheres taken 
from three orthogonal directions have been shown in an 
earlier study 

To investigate emission properties and resonance life- 
times of the deformed microspheres, we use frustrated 
total internal reflection at the surface of a prism with 
refractive index n = 1.7 to launch individual traveling- 
wave WG modes near the x-y plane of the spheroid (Fig. 
^ inset) 0. The initial angle of incidence xo with which 
the WG modes are launched in the microsphere can be 
controlled by adjusting the angle of incidence ip in the 
prism. 

WG mode spectra were obtained by measuring the 
far-field emission intensity as a function of the excita- 
tion wavelength. Since the stem holding the resonator 
is a strong leakage pathway, the observed long-lived WG 
modes must have internal field patterns that do not over- 
lap with the stem region. This is made possible by the 
fact that the resonator is slightly flattened in the z di- 
rection (along the stem axis). Experiments as well as 
ray simulations show that this deformation stabilizes the 
rays in the vicinity of the x-y plane. Our direction- 
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FIG. 1: (a)-(d): Far-field emission patterns of WG modes. 
Insets: bottom view of the resonators showing the progression 
of shapes in the x-y cross-section, from which we determine 
e = 6.7% (a), e = 2.2% (b), and e = 1.2% (c,d). WG modes 
were launched at sinxo « 1 in (a)-(c), and at sinxo ~ 0.8 
in (d); x is the internal angle of incidence with respect to 
the surface normal, (e)-(h): the spectra corresponding to the 
modes in (a)-(d), from which we deduce the Q-factors. Inset 
to (d): stable (solid) and unstable (dashed) 4-bounce orbits 
in a quadrupole at e = 1.2%. 



ality measurements correspondingly were performed in 
this plane. By scanning the detector while keeping the 
resonator fixed, we recorded WG mode spectra at var- 
ious angles, (j) (from the x-axis), to construct the far- 
field emission pattern at a given microsphere deforma- 
tion e. The microsphere was then reheated to reduce e 
and the far-field emission pattern was measured again, 
leading to the evolution shown in Figure ^ The exci- 
tation beam is s-polarized, and the detection scheme is 
polarization-insensitive. The far-field emission patterns 
are independent of the polarization of the excitation laser 
beam. Spheroids with similar e exhibit qualitatively the 
same behavior in terms of mode spectra, Q-factors, and 
emission patterns. 

Beginning with a strongly deformed microsphere, the 
emission was observed to have a strong peak at = 45° 
in the far field [see Fig. ^(a)]. Since clockwise traveling 
waves were excited, as indicated by the white arrow in 



the sphere image in Fig. ^(8-), this far-field emission di- 
rection corresponds to light escaping tangential from the 
region at = —45° on the surface of the microsphere, 
shown by the solid arrow in the sphere image. A much 
smaller peak (5% of the large peak height) in the far- 
field emission pattern was observed at cj) ~ 135°, corre- 
sponding to light escaping from the region at ^ = 45° on 
the microsphere surface (dotted arrow) . As the deforma- 
tion was reduced, the far-field emission peak at (/) = 135° 
grew to about one quarter of the 4> = 45°-peak height in 
Fig. n(b). In Fig. ^(c), the two peaks reached nearly 
equal intensity. The bright emission spots at = ±45° 
on the resonator surface can be viewed directly with a 
CCD camera. For reference, we call the pattern in Fig. 
n (a) asymmetric and the pattern in Fig. ^ (c) symmet- 
ric (around = 90°). Figures (e)-(g) also show that as 
the deformation was reduced, the Q-factor of the relevant 
WG modes increased by nearly four orders of magnitude. 

The measurements discussed thus far were performed 
with input condition xo ~ 90° (sinxo ~ !)■ Figure 
(d) shows the far-field emission pattern from the same 
microsphere as Fig. ^ (c) but with light launched at 
sinxo ~ 0.8. The emission pattern is nearly identical to 
Fign(c) where sinxo ~ 1, although the corresponding 
Q-factor is nearly two orders of magnitude smaller [Figs. 
[11(g) and (h)]. 

The observed emission patterns in Fig. n(a)-(d) dis- 
agree with the intuitive expectation that WG modes in 
an oval resonator should preferentially emit tangential to 
the points of highest curvature, into the far-field direc- 
tion <j) = ±90° . This is what one obtains when modeling 
our spheroids as triaxial ellipsoids, whose internal ray 
dynamics exhibits no chaos j^l , independently of the axis 
ratios. 

The observed emission pattern of the most deformed 
microsphere. Fig. ^(a), is well explained by a ray model. 
The peak at <^ = 45° can be attributed to an effect 
known as dynamical eclipsing 0|. In the ray model, a 
WG mode corresponds to light rays trapped close to the 
perimeter of the dielectric resonator by total internal re- 
flection, which prevents light escape unless a critical angle 
Xc = arcsin 1/n = 43.6° is reached, where n = 1.45 is the 
refractive index of the fused silica. At small but finite 
e, any oval can be approximated by the first terms of 
a multipole expansion, which after proper choice of ori- 
gin is quadrupolar. In this limit, a stable 4-bounce orbit 
shaped like a diamond forms with its sharp vertices at the 
highest-curvature points of the resonator with an angle of 
incidence Xi ~ 45° [inset to Fig. Efd)]. This creates "is- 
lands of stability" in phase space rays launched near 
the diamond orbit will retain similar reflection points and 
angles even if they do not close onto themselves after four 
bounces. At the deformation used in Fig. n](a), most of 
the phase space supporting the WG mode in question 
is chaotic. However, chaotic rays cannot penetrate the 
4-bounce islands; this means chaotic WG rays are pre- 
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vented from refractively escaping at the points of highest 
curvature, because X4 ~ Xc- Instead, as trajectories flow 
around the islands in one direction, refractive escape oc- 
curs near </> = —45° on the resonator surface. From this 
circulation around the islands, asymmetric far-field emis- 
sion patterns result as a hallmark of refractive escape. 
This is what we observe in Fig. ^(a). 

This ray mechanism was originally proposed and tested 
for effectively 2D systems where islands of stability 
and chaotic regions are mutually exclusive. In our 3D 
spheroids, the same diamond-shaped stable orbit exists 
in or near the x-y plane, and refractively escaping rays in 
its vicinity have lifetimes that correspond to < 10 round 
trips in the resonator, translating to a Q-factor near 10**. 
During this time, a given ray can be considered as mov- 
ing in a cross-sectional plane that may be inclined and 
slowly precessing around the x-axis. The variation of 
cross-sectional shape, Ae, probed by such rays causes no 
significant differences in the size of the 4-bounce islands 
[iQj . Therefore, rays in the 3D resonator behave as in 
a corresponding fixed 2D resonator during a relatively 
short time preceding an escape event, and in particular 
show dynamical eclipsing. 

However, refractive escape ceases to be the dominant 
leakage mechanism at the very small deformation used 
in Figs, m (c) and (d), since in this case chaotic diffu- 
sion is no longer effective and the stable diamond orbit 
itself also becomes fully confined by total internal reflec- 
tion. Those classical rays that do escape refractively, 
moreover, are not able to produce the symmetric emis- 
sion patterns observed in Figs, ^(c) and (d) 0|. Hence, 
the emission process in these cases is governed by evanes- 
cent escape, i.e. a tunnel coupling between internal and 
radiation fields. 

Tunneling is the only decay mechanism in an ideal 
sphere where each WG mode in the plane of excitation 
corresponds to a unique azimuthal angular momentum 
number m with respect to the z-axis, which is semiclas- 
sically related to xo by sinxo — m/(nkR), using the fact 
that the modes of interest remain close to the x-y plane 
and hence have total angular momentum I k, m. At 
kR ^ I, the tunneling escape rate is negligible when 
sinxo ~ 1 but accelerates exponentially as sinxo ap- 
proaches sinxc- In a deformed microsphere, the an- 
gle of incidence x varies as a well-defined function of 
provided it covers an angular momentum range where 
chaos is neglibible [l^ . A circulating ray with varying 
x{4') > Xc will then escape with exponentially strong se- 
lectivity near the minima of x(0). 

The emission locations and directions observed in the 
experiment indicate that the minima of xi4') near the 
corners of the unstable rectangular 4-bounce orbit [inset 
to Fig. n(d)]. We determined from ray tracing that this 
only occurs near but above Xc- Thus, the symmetric far- 
field patterns in Figs, ^(c) and (d) are a tunneling probe 
for the minima of the ray-optical x{4') in the vicinity of 



the critical angle Xc- To further corroborate this inter- 
pretation, recall that the emission patterns as shown in 
Figs, n (c) and (d) are insensitive to the initial angle 
of incidence, xoi as they should be if the detected light 
originates near sinx ~ sinxc- The question remains how 
the light coupled into the resonator is able to reduce its 
angular momentum from a high value corresponding to 
sinxo to the much lower sinXc- Our observation shows 
that such a dynamical process is present in the spheroid, 
but does not unambiguously reveal its mechanism. A 
possible mechanism is Arnol'd diffusion, a phase-space 
transport process that crucially depends on the fact that 
the resonators are in fact three-dimensional 0, with no 
axial symmetry. Even without resolving this question, 
we can nevertheless conclude that the phase space re- 
gion near sinXc is responsible for the observed emission 
patterns. 

Our initial assumption of a 2D model for the emission 
process, justified by the approximate 2D nature of the 
3D ray trajectories for intermediate times as discussed 
earlier, is therefore consistent with the fact that no long- 
time properties of the internal ray dynamics are needed 
to explain the emission directions. With the same justi- 
fication, we also performed numerical wave calculations 
by modeling the shape felt by modes in the x-y plane 
as the cross section of a quadrupolar cylinder. Figures El 
(a) and (b) plot the intensity patterns of two WG modes 
under traveling- wave excitation 0, 0| . The calculation 
assumes s-polarized light and uses size parameters near 
kR^ 113 Q. 

The emission patterns in Figs. 121(a) and (b) agree qual- 
itatively with the asymmetric and the symmetric emis- 
sion patterns shown in Figs. n](b) and (d), respectively. 
To understand the physical mechanism for the emission 
in Fig. 121 (b), we plot in Fig. [21 (c) the distribution 
of angular momentum numbers for the WG mode in (b). 
Figure|21(c) shows negligible overlap with the window for 
refractive escape below the critical angle Xc (correspond- 
ing to a critical angular momentum nic ~ kR « 113). 
This indicates that the emission from this mode is due 
to tunneling escape. By comparison, both refractive and 
tunnehng escape contribute to the emission in Fig. 121(a), 
with the refractive escape playing a dominant role. 

The peak splitting of width Am « 4 in Fig. (21 (c) is 
a straightforward consequence of the oval deformation: 
as the wave circulates around the resonator with varying 
radius between Rmin and Rmax, its angular momentum 
oscillates but has high probabitity of being near its ex- 
trema, given by the extrema of m f» r((/<)nfci? sin x(0) 
over the surface. The fact that the minima of x{4') do 
not occur a,t cj) — 0° , 180° on the surface is due to the 4- 
bounce island structure in the vicinity of sinx = sin 45°. 

For comparison, Fig. (21 (c) also shows the distribu- 
tion of angular momentum number of an ultra-high Q 
WG mode with emission directionality as expected for 
the ellipse [see Fig. 121(d)]. The mode is confined at a 
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FIG. 2: Top: Intensity pattern of WG modes with parameters 
kR = 112.063, £ = 6.5% (a) and kR = 112.452, e = 3.4% (b). 
The Q-factors are 4 x 10* (a) and 2 x 10^ (b), and the overaU 
directionality agrees with that in Fig. (b) and (d) , respec- 
tively. The values of e differ from Fig. Q because the simu- 
lation approximates the unknown details of the experimental 
shape by a pure quadrupole 0|. Bottom: solid histogram 
in (c) shows the angular-momentum distribution of the mode 
in (b). A very high-Q WG mode at the same deformation 
(c, dashed), with kR = 113.229 and Q = 3 x lO" radiates 
predominantly into tp = 90° , as shown in the far-field pattern 
(d). 



high m, corresponding to sinxo ~ 0.91. Note that here 
the tunneUng loss is neghgible compared with scattering 
or absorption loss of the material that limits the actual 
Q-factor of fused silica microspheres to of order 10^ . 
The two modes in Fig. [21(c) have practically no overlap in 
angular momentum space. The qualitative difference in 
the far-field patterns of these two modes further confirms 
that the mode shown in Fig. ^ (c), which was excited at 
sinxo ~ 1, is not confined to sinxo ~ 1- 

In essence, our experiment exploits the peculiar fea- 
ture that for fused silica the critical angle lies near the 
phase-space islands corresponding to the stable 4-bounce 
orbit. These islands are the dominant structure in the 
WG region and persist even at small e where most of the 
phase space is filled with non-chaotic, "regular" trajecto- 
ries. The island formation is a non-perturbative conse- 
quence of the breakdown of conservation laws |^, which 
qualitatively distinguishes the generic quadrupolar shape 
from an ideal ellipse even though the two shapes dif- 
fer only to second order in e. The emission character- 
istics of a silica spheroid in the seemingly trivial small-e 
regime will continue to be strongly affected by ray pat- 
terns that wrap around the perimeter in approximately 



four bounces. 

Whether any non-perturbative structure in the ray dy- 
namics can be resolved by the wave field, depends on kR 
[lfij |: a directionality measurement will be able to distin- 
guish the peaks a,t cj) — 90° ±45° if the conjugate angular 
momentum m satisfies the uncertainty relation Am > 
1/2 with Acj) « 7r/4. This implies Am > 2/7r w 1, which 
in our spheroids translates to a fluctuating angle of in- 
cidence of Asinx — Am/{nkR) « 10"'^. Even the least 
deformed resonator (e « 1%) exceeded this estimated 
resolution threshold significantly. 

It is truly remarkable that this intricate interplay be- 
tween ray and wave dynamics for which chaos plays no 
dominant role, can be exploited to engineer WG modes 
that can feature both high-Q and highly directional emis- 
sion. This makes WG resonators with small deformation 
highly promising for a variety of applications |l7t Il8l |. 
such as microlasers, nonlinear optics, and quantum in- 
formation processing. 

This work was supported in part by NSF under grants 
No. DMR9733230 and No. DMR0201784. 
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